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^ ' Abstract 

I— I. Récent results about the coexistence of ferromagnetism and unconventional su- 

Q , perconductivity with spin-triplet Cooper pairing are reviewed on the basis of the 

^ I quasi-phenomenological Ginzburg-Landau theory. The superconductivity in the 
mixed phase of coexistence of ferromagnetism and unconventional superconductivity 

^ ' is triggered by the spontaneous magnetization. The mixed phase is stable whereas 

^ ■ the other superconducting phases that usually exist in unconventional supercon- 

^ • ductors are either unstable or metastable at relatively low températures in a quite 

G ■ narrow domain of the phase diagram and the stability properties are determined by 

. the particular values of Landau parameters. The phase transitions from the normal 

Q , phase to the phase of coexistence is of first order while the phase transition from 

O I the ferromagnetic phase to the coexistence phase can be either of first or second 

I order depending on the concrète substance. The Cooper pair and crystal anisotropy 

^ I are relevant to a more précise outline of the phase diagram shape and reduce the 

' degeneration of the ground states of the system. The results are discussed in view 

^ ' of application to itinérant ferromagnetic compounds as UGe2, ZrZn2, URhGe. 
\0 ■ 

O ■ 

m : 

o . 1 Introduction 

■ 

a : 

V . 1.1 Notes about unconventional superconductivity 
'fi ' 

O ! The phenomenon of unconventional Cooper pairing of fermions, i.e., the formation of 
! Cooper pairs with nonzero angular momentum was theoretically predicted ilj in 1959 
.!-H ! as a mechanism of superfluidity in Fermi liquids. In 1972 the same phenomenon - un- 
^ ! conventional superfluidity due to a p-wave (spin triplet) Cooper pairing of ^He atoms, 
! was experimentally discovered in the mK range of températures; for détails and the- 
oretical description see Refs. [21 El El- Note that, in contrast to the standard s-wave 
pairing in usual (conventional) superconductors where the électron pairs are formed by 
an attractive electron-electron interaction due to a virtual phonon exchange, the widely 
accepted mechanism of the Cooper pairing in superfluid ^He is based on an attractive 
interaction between the fermions (^He atoms) as a resuit of a virtual exchange of spin 
fluctuations. Certain spin fluctuation mechanisms of unconventional Cooper pairing of 
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électrons are assumed also for the discovered in 1979 heavy fermion superconductors (see, 
e.g., Refs. iniElIZI) as well as for some classes of higli-temperature superconductors (see, 

e.g., Refs. jHiiniiiniiniiiaiiniiiiiiisiiin]). 

The possible superconducting phases in unconventional superconductors are described in 
the framework of the gênerai Ginzburg-Landau (GL) effective free energy functional JH] 
with the help of the symmetry groups theory. Thus a variety of possible superconduct- 
ing orderings were predicted for différent crystal structures [T71 HHl HHl 1201 122] ■ A 
detailed thermodynamic analysis fH] ITH] of the homogeneous (Meissner) phases and a 
renormalization group investigation JT] of the superconducting phase transition up to 
the two-loop approximation have been also performed (for a three-loop renormalization 
group analysis, see Réf. [22]; for effects of magnetic fluctuations and disorder, see |2Hl25]). 
We shall essentially use thèse results in our présent considération. 

1.2 Expérimental prédictions 

In 2000, experiments j2n| at low températures (T ~ 1 K) and high pressure (P ~ 1 
GPa) demonstrated the existence of spin triplet superconducting states in the metallic 
compound UGe2. The superconductivity is triggered by the spontaneous magnetization 
of the ferromagnetic phase that occurs at much higher températures. It coexists with 
the superconducting phase in the whole domain of its existence below T ~ 1 K; see also 
experiments from Refs. [2ZII2H]! and the discussion in Réf. [201 • ^he same phenomenon 
of existence of superconductivity at low températures and high pressure in the domain 
of the (T, P) phase diagram where the ferromagnetic order is présent was observed in 
other ferromagnetic metallic compounds (ZrZn2 IHO] and URhGe [SI]) soon after the 
discovery ^Hl of superconductivity in UGe2. 

In superconducting ternary and Chevrel compounds the influence of magnetic order on 
superconductivity is also substantial (see, e.g., [321 E31 EU ESI) t>ut in the newly found 
ferromagnetic substances the phase transition température (Tf) to the ferromagnetic state 
is much higher than the phase transition température (Tps) from ferromagnetic to a mixed 
state of coexistence of ferromagnetism and superconductivity. For example, in UGe2, 
Tps = 0.8 K while the critical température of the phase transition from paramagnetic to 
ferromagnetic state in the same material is = 35 K [201 12Z|- It can be assumed that 
in thèse substances the material parameter Tg defined as the usual critical température 
of the second order phase transition from normal to uniform (Meissner) supercondicting 
state in a zéro external magnetic field is much lower than the phase transition température 
Tps- The above mentioned experiments on the compounds UGe2, URhGe, and ZrZn2 do 
not give any évidence for the existence of a standard normal-to-superconducting phase 
transition in a zéro external magnetic field. 

It seems that the superconductivity in the metallic compounds mentioned above always 
coexists with the ferromagnetic order and is enhanced by it. In thèse Systems, as claimed 
in Réf. [20], the superconductivity probably arises from the same électrons that create 
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the band magnetism and can be most naturally understood rather as a triplet than spin- 
singlet pairing phenomenon. Metallic compounds UGe2, URhGe, and ZrZn2, are itinérant 
ferromagnets. An unconventional superconductivity is also suggested jSEl as a possible 
outcome of récent experiments in Fe j^, in which a superconducting phase has been 
discovered at températures below 2 K and pressures between 15 and 30 GPa. There 
both vortex and Meissner superconductivity phases [ÏÏTj are found in the high-pressure 
crystal modification of Fe with a hexagonal close-packed lattice for which the strong 
ferromagnetism of the usual bec iron crystal probably disappears (SEI- It can be hardly 
claimed that in hexagonal Fe the ferromagnetism and superconductivity coexist but the 
clear évidence for a superconductivity is also a remarkable achievement. 

The reasonable question whether thèse examples of superconductivity and coexistence 
of superconductivity and ferromagnetism are bulk or surface effects can be stated. The 
earlier experiments performed before 2004 do not answer this question. Récent experi- 
ments jHH] show that surface superconductivity appears in ZrZn2 and its présence dépends 
essentially on the way of préparation of the sample. But in our study it is important that 
bulk superconductivity can be considered well established in this substance. 

1.3 Ferromagnetism versus superconductivity 

The important point in ail discussions of the interplay of superconductivity and ferro- 
magnetism is that a small amount of magnetic impurities can destroy superconductivity 
in conventional (s-wave) superconductors by breaking up the (s-wave) électron pairs with 
opposite spins (paramagnetic impurity effect [^). In this aspect the phenomenological 
arguments [10] and the conclusions on the basis of the microscopic theory of magnetic 
impurities in s-wave superconductors [HH! are in a complète agreement with each other; 
see, e.g., Refs. [221 ESI EU ES] • In fact, a total suppression of conventional (s-wave) su- 
perconductivity should occur in the présence of an uniform spontaneous magnetization 
M, i.e., in a standard ferromagnetic phase [10]. The physical reason for this suppression 
is the same as in the case of magnetic impurities, namely, the opposite électron spins in 
the s-wave Cooper pair turn over along the vector M in order to lower their Zeeman 
energy and, hence, the pairs break down. Therefore, the ferromagnetic order can hardly 
coexist with conventional superconducting states. Especially, this is valid for the coex- 
istence of uniform superconducting and ferromagnetic states where the superconducting 
order parameter tp{x) and the magnetization M do not dépend on the spatial vector x. 

But yet a coexistence of s-wave superconductivity and ferromagnetism may appear in un- 
common materials and under quite spécial circumstances. Furthermore, let us emphasize 
that the conditions for the coexistence of nonuniform ( "vertex" , "spiral" , "spin-sinosoidal" 
or "helical" 32^ ,33j) superconducting and ferromagnetic states are less restrictive than 
those for the coexistence of uniform superconducting and ferromagnetic orders. Coexis- 
tence of nonuniform phases has been discussed in détails, both experiment and theory, in 
ternary and Chevrel-phase compounds where such a coexistence seems quite likely; for a 
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comprehensive review, see, for example, Refs. [H^ IHHl ^^ IHïïl HT]. 

In fact the only two superconducting Systems for which the expérimental data allow 
assumptions in a favor of a coexistence of superconductivity and ferromagnetism are the 
rare earth ternary boride compound ErRli4B4 and the Chervel phase compound HoMoeSg; 
for a more extended review, see Refs. In thèse compounds the phase of coexistence 

appears in a very narrow température région just below the Curie température Tf of the 
ferromagnetic phase transition. At lower températures the magnetic moments of the rare 
earth 4/ électrons become better aligned, the magnetization increases and the s-wave 
superconductivity pairs formed by the conduction électrons disintegrate. 

1.4 Unconventional superconductivity triggered by ferromag- 
netic order 

We shall not consider ail important aspects of the long standing problem of coexistence 
of superconductivity and ferromagnetism rather we shall concentrate our attention on the 
description of the newly discovered coexistence of ferromagnetism and unconventional 
(spin-triplet) superconductivity in the itinérant ferromagnets UGe2, ZrZn2, and URhGe. 
Here we wish to emphasize that the main object of our discussion is the superconductivity 
of thèse compounds and at a second place in the rate of importance we put the problem of 
coexistence. The reason is that the existence of superconductivity in itinérant ferromag- 
nets is a highly nontrivial phenomenon. As noted in Réf. |32] the superconductivity in 
thèse materials appears to be difficult to explain in terms of previous théories [321 IHSI 
and requires new concepts to interpret the expérimental data. 

We have already mentioned that in ternary compounds the ferromagtetism comes from 
the localized 4/ électrons while the s-wave Cooper pairs are formed by conduction élec- 
trons. In UGe2 and URhGe the 5/ électrons of U atoms form both superconductivity 
and ferromagnetic order |26] lÏÏT] . In ZrZn2 the same double rôle is played by the Ad élec- 
trons of Zr. Therefore, the task is to describe this behavior of the band électrons at a 
microscopic level. One may speculate about a spin-fluctuation mediated unconventional 
Cooper pairing as is in case of ^He and heavy fermion superconductors. Thèse important 
issues have not yet a reliable answer and for this reason we shall confine our considération 
to a phenomenological level. 

In fact, a number of reliable expérimental data as the cohérence length and the supercon- 
ducting gap measurements jSniEZUHHIHni are in favor of the conclusion about a spin-triplet 
Cooper pairing in thèse metallic compounds, although the mechanism of pairing remains 
unclear. We shall essentially use this reliable conclusion. This point of view is consistent 
with the expérimental observation of coexistence of superconductivity only in a low tem- 
pérature part of the ferromagnetic domain of the phase diagram (T, P), which means that 
a pure (non-ferromagnetic) superconducting phase is not observed, a circumstance, that 
is also in favor of the assumption of a spin-triplet superconductivity. Our investigation 
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leads to results which confirm this gênerai picture. 

On the basis of the expérimental data and conclusions presented for the first time in 
Refs. jSni lin] and shortly afterwards confirmed in Refs. j23 EH IHOl IHIj one may reliably 
accept that the superconductivity in thèse magnetic compounds is considerably enhanced 
by the ferromagnetic order parameter M and, perhaps, it could not exist without this 
"mechanism of ferromagnetic trigger," or, in short, "iVf-trigger" ; see Refs. |^ I45| Hïïj 
where this concept has been introduced for the first time. The trigger phenomenon is 
possible for spin-triplet Cooper pairs where the électron spins point parallel to each other 
and their turn along the vector of the spontaneous magnetization M does not produce 
a breakdown of the spin-triplet Cooper pairs but rather stabilizes them and, perhaps, 
stimulâtes their création. We shall describe this phenomenon at a phenomenological 
level. 

1.5 Phenomenological studies 

A phenomenological theory that explains the coexistence of ferromagnetism and uncon- 
ventional spin-triplet superconductivity of GL type has been developed recently in [intITTj 
where possible low-order couplings between the superconducting and ferromagnetic order 
parameters are derived with the help of gênerai symmetry group arguments. On this 
basis several important features of the superconducting vortex state of unconventional 
ferromagnetic superconductors were demonstrated 1^ . 

In our review we shall follow the approach from Refs. |13 El to investigate the condi- 
tions for the occurrence of the Meissner phase and to demonstrate that the présence of 
ferromagnetic order enhances the p-wave superconductivity. We also establish the phase 
diagram of ferromagnetic superconductors in a zéro external magnetic field and show that 
the phase transition to the superconducting state can be either of first or second order 
depending on the particular substance. We confirm the prédictions made in Refs. |1SII1Z| 
about the symmetry of the ordered phases. 

In our study we use the mean-field approximation jî^j and known results about the 
possible phases in nonmagnetic superconductors with triplet (p-wave) pairing [THl HH 1121 
10]. Our results [SllinillOi show that taking into account the anisotropy of the spin-triplet 
Cooper pairs modifies but does not drastically change the thermodynamic properties of 
the coexistence phase, especially in the température domain above the superconducting 
critical température Tg. The effect of crystal anisotropy is similar but we shall not make 
an overall thermodynamic analysis of this problem because we have to consider concrète 
Systems and crystal structures for which there is no enough information from 

experiment to make conclusions about the parameters of the theory. Our results confirm 
the gênerai concept that the anisotropy reduces the degree of ground state degeneration, 
and depending on the symmetry of the crystal, picks up a crystal direction for the ordering. 

There exists a formai similarity between the phase diagram we obtain and the phase di- 
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agram of certain improper ferroelectrics (HUl IHH 1^21 IMl IMl 1^^] • We shall make use of 
the concept in the theory of improper ferroelectrics, where the trigger of the primary 
order parameter by a secondary order parameter (the electric polarization) has been ini- 
tially introduced and exploited; see Réf. [SnilMlISSl- The mechanism of the M-triggered 
superconductivity in itinérant ferromagnets |S1 SU HE] is formally identical to the mech- 
anism of appearance of structural order triggered by the electric polarization in improper 
ferroelectrics (see, e.g., Refs. jHni EU ESI)- 

Our investigation is based on the GL free energy functional of unconventional ferromag- 
netic superconductors presented in Sec. 2.1 and we shall establish the uniform phases 
which are described by it. More information about the justification of this investigation 
is presented in Sec. 2.2. We work with a quite gênerai GL free energy and the problem is 
that there is no enough information about the values of the parameters of the model for 
concrète compounds (UGe2, URhGe, ZrZn2) where the ferromagnetic superconductivity 
has been discovered. On the one hand the lack of information makes impossible a de- 
tailed comparison of the theory to the available expérimental data but on the other hand 
our results are not bound to one or more concrète substances and can be applied to any 
unconventional ferromagnetic superconductor. In Sec. 3 we discuss the phases in non- 
magnetic unconventional superconductors. In Sec. 4 the M-trigger effect |^ EE] will 
be described when only a linear coupling of the magnetization M to the superconducting 
order parameter ip is considered; here the spatial dependence of order parameters and ail 
anisotropy effects are ignored. In Sec. 5 we analyze the influence of quadratic coupling of 
magnetization to the superconducting order parameter on the thermodynamics of the fer- 
romagnetic superconductors. The application of our results to expérimental (T, P) phase 
diagrams is discussed in Sec. 5.3. In Sec. 6 the anisotropy effects are outlined. Our main 
attention is focussed on the Cooper-pair anisotropy. Note, that certain types of crystal 
anisotropy may produce more than one ferromagnetic phase |56t \57\ l58] but here we shall 
not dwell on this interesting topic. In Sec. 7 we summarize our conclusions. 

2 Ginzburg-Landau free energy 

Following Refs. jlHl IIHI HH HHl HEl HZ] in this Chapter we discuss the phenomenologi- 
cal theory of spin-triplet ferromagnetic superconductors and justify our considération in 
Sections 3-6. 



2.1 Model 



The gênerai GL free energy functional, we shall use in our analysis, is 
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where the free energy density f{ip,M) ( hereafter called "free energy") of a spin-triplet 
ferromagnetic superconductor is a sum of five terms 113 ^] , namely, 

/(^, M) = fsiij) + fUM) + fiiij, M) + ^- B.M. (2) 

In Eq. (2) the three dimensional complex vector ip = {ipj'jj = 1,2,3} represents the 
superconducting order parameter, B = {H + AtcM) = V x A is the magnetic induction; 
H is the external magnetic field, A = {Aj;j = 1,2,3} is the magnetic vector potentiaL 
The last two terms on r.h.s. of Eq. (2) are related with the magnetic energy which includes 
both diamagnetic and paramagnetic effects in the superconductor; see, e.g., pSi |55]. 

The energy part fsii') in Eq. (2) describes the superconductivity for H = M = 0. It can 
be written in the form 

Here 

where a summation over the indices i,j = 1, 2, 3 is assumed and the symbol 

n . = -ih^ + '^A. (5) 

OXi c 

dénotes a covariant differentiation. In Eq. (3), > and = as{T — T,), where is 
a positive material parameter and Tg is the critical température of the standard second 
order phase transition which may occur at H = A4 = 0; H = \H\, and Ai = \M\. The 
quantities Ug and Vg describe the anisotropy of the spin-triplet Cooper pair and the crystal 
anisotropy, respectively, [IHIIII]. Parameters Kj, (j = 1,2,3) in Eq. (4) are related with 
the effective mass tensor of anisotropic Cooper pairs [TH] . 

The superconducting part (3) of the free energy f{ip, M) is derived from symmetry group 
arguments and is independent of particular microscopic models; see, e.g., Refs. [T^ IH]. 
According to classifications [IHl Ej the p-wave superconductivity in the cubic point group 
Oh can be realized through one-, two-, and three-dimensional représentations of the order 
parameter. The expressions (3) and (5) incorporate ail three possible cases. The coef- 
ficients bg, Ug, and Vg in Eq. (3) are différent for weak and strong spin-orbit couplings 
but in our investigation they are considered as undetermined material parameters which 
dépend on the particular substance. 

The free energy of a standard isotropic ferromagnet is given by the term fp{M) in Eq. (2), 
/f(M) =Cff2 I V.^.f + «/(T;)M2 + |M^ (6) 
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where Vj = d/dxj and hj > 0. The quantity af{T'j) = «/(T — Tj.) is expressed by the 
material parameter aj > and the température Tj which is différent from the critical 
température Ty of the ferromagnet and this point will be discussed below. We have 
already added a négative term (— 27rA^^) to the total free energy /(?/', Ai") and that is 
obvions by setting H = in Eq. (2). The négative energy (— 27rA1^) should be added to 
f^{M). In this way one obtains the total free energy fp{M) of the ferromagnet in a zéro 
external magnetic field that is given by a modification of Eq. (6) according to the rule 

Maf) = f^[ajiT'j)^ajiTj)], (7) 

where a/ = a/(T — Tf) and 

r, = T; + | (8) 

is the critical température of a standard ferromagnetic phase transition of second order. 
This scheme was used in studies of rare earth ternary compounds [221 EHl EOl El] • Alterna- 
tively [021 5 from the beginning the total ferromagnetic free energy fpiaf, M) 

as given by Eqs. (6) - (8) but in this case the magnetic energy included in the last two 
terms on r.h.s. of Eq. (2) should be replaced with H'^/Sn. Both approaches are équivalent. 

The interaction between the ferromagnetic order parameter M and the superconducting 
order parameter ip is given by jlSl E2| 

fi{^, M) = t^oM.{^ X r) + ÔM^\^\'. (9) 

The 7o-term in the above expression is the most substantial for the description of ex- 
perimentally found ferromagnetic superconductors |17j and the ôM'^\tjj\'^-teTm makes the 
model more realistic in the strong coupling limit as it gives the opportunity to enlarge the 
phase diagram including both positive and négative values of the parameter a^. In this 
way the domain of the stable ferromagnetic order is extended down to zéro températures 
for a wide range of values of material parameters and the pressure P, a situation that 
corresponds to the experiments in ferromagnetic superconductors. 

In Eq. (9) the coupling constant 7o > can be represented in the form 70 = AttJ, where 
J > is the ferromagnetic exchange parameter |17]. In gênerai, the parameter ô for 
ferromagnetic superconductors may take both positive and négative values. The values 
of the material parameters (T^, Tf, as, af, hg, Ug, Vg, bf, Kj, 70 and ô) dépend on the 
choice of the concrète substance and on thermodynamic parameters as température T 
and pressure P. 



2.2 Way of treatment 

The total free energy (2) is a quite complex object of theoretical investigation. The 
possible vortex and uniform phases of the model cannot be investigated within a single 
calculation, rather one should focus on concrète problems. In Réf. jUj the vortex phase 
was discussed with the help of the criterion [HBl for a stability of this state near the phase 
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transition line Tc2{H); see also, Réf. [03]. In case oî H = one should apply the same 
criterion with respect to the magnetization Ai for small values of near the phase 
transition line Tc2{Ai) as performed in Réf. 

Here we shall be interested in the uniform phases when the order parameters ip and M 
do not dépend on the spatial vector x E V, where V is the volume of the superconductor. 
We shall restrict our analysis to the considération of the coexistence of uniform (Meiss- 
ner) phases and ferromagnetic order. We shall make a detailed investigation in order to 
show that the main properties of the uniform phases can be well determined within an 
approximation when the crystal anisotropy is neglected. Even, some of the main features 
of the uniform phases in unconventional ferromagnetic superconductors can be reliably 
outlined when the Cooper pair anisotropy is neglected, too. 

The assumption of a uniform magnetization M is always reliable outside a quite close 
vicinity of the magnetic phase transition and under the condition that the superconduct- 
ing order parameter ip is also uniform, i.e. that vortex phases are not présent at the 
respective température domain. This conditions are directly satisfied in type I supercon- 
ductors but in type II superconductors the température should be sufïiciently low and 
the external magnetic field should be zéro. Nevertheless, the mentioned conditions for 
type II superconductors may turn insufficient for the appearance of uniform supercon- 
ducting states in materials with quite high values of the spontaneous magnetization. In 
such situation the uniform (Meissner) superconductivity and, hence, the coexistence of 
this superconductivity with uniform ferromagnetic order may not appear even at zéro 
température. Up to now type I unconventional ferromagnetic superconductors have not 
been found whereas the expérimental data for the recently discovered compounds UGe2, 
URhGe, and ZrZn2 are not enough to conclude definitely either about the lack or the 
existence of uniform superconducting states at low and ultra-low températures. 

If real materials can be modelled by the gênerai GL free energy (1) - (9), the ground 
state properties will be described by uniform states, which we shall investigate. The 
problem about the availability of such states in real materials at finite températures is 
quite subtle at the présent stage of research when the expérimental data are not enough. 
Recently in Réf. [0^] an expérimental évidence was given for the coexistence of uniform 
superconductivity and ferromagnetism in UGe2. Thus we shall assume that uniform 
phases may exist in some unconventional ferromagnetic superconductors. Moreover, we 
have to emphasize that thèse phases appear as solutions of the GL équations corresponding 
to the free energy (1) - (9). Thèse arguments completely justify our study. 

In case of a strong easy axis type of magnetic anisotropy, as is in UGe2 [211; the overall 
complexity of mean-field analysis of the free energy f{ip, M) can be avoided by performing 
an "Ising-like" description: M = (0, 0, Ai), where Ai = ±|-/Vf| is the magnetization along 
the "2;-axis." Because of the équivalence of the "up" and "down" physical states (±-/W) 
the thermodynamic analysis can be performed within the "gauge" > 0. When the 
magnetic order has a continuons symmetry we can take advantage of the symmetry of 
the total free energy /{ip, M) and avoid the considération of équivalent thermodynamic 
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States that occur as a resuit of the respective symmetry breaking at the phase transition 
point but have no effect on thermodynamics of the System. In the isotropic System one 
may again choose a gauge, in which the magnetization vector has the same direction as 
2;-axis {\M\ — = Ai) and this wiU net influence the generality of thermodynamic 
analysis. Here we shall prêter the alternative description within which the ferromagnetic 
state may appear through two équivalent up and down domains with magnetizations M. 
and {—Ai), respectively. 

Wc shall make the mean-field analysis of the uniform phases and the possible phase 
transitions bctwcen such phases in a zéro external magnetic field (H = 0), when the 
crystal anisotropy is neglcctcd {vs = 0). The only exception will be the considération in 
Sec. 3, where we briefly discuss the nonmagnetic superconductors (Ai ^ 0). 

We shall use notations in which the number of parameters is reduced. For this reason we 
introduce 

b={bs+Us + V,) (10) 

and redefine the order parameters and the other quantifies in the foUowing way: 

= 6V4^^. = ^ M = b^^M , (11) 
as ^ af Us Vs 

^=17' ^ = T' 
Vb y/b} b b 

7o S 
7 = 7— T7Ï' 71 



61/2^,1/4' - (66^)1/2- 

Having in mind that the order parameters ip and M are considered uniform and using 
Eqs. (10) and (11), we can write the free energy density f{ip,M) = F{ip, M)/V in the 
form 

f{iP,M) = r02 + ^(/)^ + 270i(/)2Msin(^2-^i) (12) 

-2w[(Pl4sm''{e2 - ei) + (i)l(Plsm\ei - es) 

+02</,2sin2(^2 - ^3)] 

Note, that in the above expression the order parameters tjj and M are defined per unit 
volume. 

The equilibrium phases are obtained from the équations of state 

df(l^o) 



, (13) 



where the séries of symbols /i can be defined as, for example, // = {na} = {M, 01, 03, 
^1, ...,^3); /iq dénotes an equilibrium phase. The stability matrix F of the phases /iq is 
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defined by 

Fi,.) ^ iFM) - |5g . (14) 

An alternative treatment can be do ne in terms of real {ipj) and imaginary [tpj) parts of 
the complex numbers ipj = V'j + ^"^J- The calculation with moduli (pj and phase angles 
6j of ipj is more simple but in cases of strongly degenerate phases some of the angles dj 
remain unspecified. Then an alternative analysis with the help of the components ipj and 
ipj should be done. 

The thermodynamic stability of the phases that are solutions of Eqs. (13) is checked with 
the help of the matrix (14). An additional stability analysis is done by the comparison 
of free énergies of phases that satisfy (13) and render the stability matrix (14) positive in 
one and the same domain of parameters {r, t, 7, 71, w, v}. This step is important because 
the complicated form of the free energy générâtes a great number of solutions of Eqs. (13) 
and we have to sift out the stable from metastable phases that correspond either to global 
or local minima of the free energy, respectively [T!^ . 

Some solutions of Eqs. (13) have a marginal stability, i.e., their stability matrix (14) is nei- 
ther positively nor negatively definite. This is often a resuit of the degeneration of phases 
with broken continuons symmetry. If the reason for the lack of a clear positive definite- 
ness of the stability matrix is precisely the mentioned degeneration of the ground state, 
one may reliably conclude that the respective phase is stable. If there is another reason, 
the analysis of the matrix (14) will be insufficient to détermine the respective stability 
property. Thèse cases are quite rare and occur for particular values of the parameters 
{r,t,7, ...}. 



3 Pure superconductivity 

Let us set M = in Eq. (12) and briefiy summarize the known results |TS1E] for the "pure 
superconducting case" when the magnetic order cannot appear and magnetic effects do not 
affect the stability of the uniform (Meissner) superconducting phases. The possible phases 
can be classified by the structure of the complex vector order parameter = [ipi, ?/'2, ^'2)- 
We shall often use the moduli vector (0i, 02,03) with magnitude = (0^ + 02 + 'PîY^'^ 
and the phase angles 9j. 

The normal phase (0,0,0) is always a solution of the Eqs. (13). It is stable for r > 0, and 
corresponds to a free energy / = 0. Under certain conditions, six ordered phases PUE] 
occur for r < 0. Here we shall not repeat the detailed description of thèse phases ITT] 
but we shall briefiy mention their structure. 

The simplest ordered phase is of type (-01,0,0) with équivalent domains: (0,-02,0) and 
(0,0,-03). Multi- domain phases of more complex structure also occur, but we shall not 
always enumerate the possible domains. For example, the "two-dimensional" phases can 
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be fully représentée! by domains of type but there are also other two types of 

domains: 0, ■ï/'s) and (0, '?/'2, V's)- As we consider the gênerai case when the crystal 
anisotropy is présent {v ^ 0), this type of phases possesses the property {ipil = {ipjl. 

The two-dimensional phases are two and have différent free énergies. To clarify this point 
let us consider, for example, the phase (t/'i, '?/'2; 0). The two complex numbers, ipi and 
ip2 can be represented either as two-component real vectors, or, equivalently, as rotating 
vectors in the complex plane. One can easily show that Eq. (12) yields two phases: a 
collinear phase, when {62 — 61) = Tik{k = 0, ±1, ...), i.e. when the vectors ipi and ip2 are 
coUinear, and another (non-collinear) phase when the same vectors are perpendicular to 
each other: (6*2 — ^1) = vr(A; + 1/2). Having in mind that |0i| = |02| = 0/\/2, the domain 
('î/'i; V'2, 0) of the collinear phase is given by (±1, 1, 0)(f)/\/2 whereas the same domain for 
the non-coUinear phase will be (±z, 1, O)0/v^. The two domains of thèse phases have 
similar représentations. 

In addition to the mentioned three ordered phases, three other ordered phases exist. For 
thèse phases ail three components ipj have nonzero equilibrium values. Two of them 
have equal to one another moduli (pj, i.e., 0i = 02 = 03- The third phase is of the 
type 01 = 02 7^ 03 and is unstable so it cannot occur in real Systems. The two three- 
dimensional phases with equal moduli of the order parameter components have différent 
phase angles and, hence, différent structure. The différence between any couple of angles 
6j is given by ±7r/3 or ±27r/3. The characteristic vectors of this phase can be of the form 
(^giTT/s^g-iTT/s^ and (e^^'^/^ e'^^'^/^ 1)0/^3. The second stable three dimensional 

phase is "real", i.e. the components ijjj lie on the real axis; {6j — Oj) = nk for any couple 
of angles 9j and the characteristic vectors are (±1, ±1, 1)0/a/3. The stability properties 
of thèse five stable ordered phases were presented in détails in Refs. ITT|. 

When the crystal anisotropy is not présent {v = 0) the picture changes. The increase 
of the level of degeneracy of the ordered states leads to an instability of some phases 
and to a lack of some noncollinear phases. Both two- and three-dimensional real phases, 
where {6j — 6j) = irk, are no more constrained by the condition 0^ = 0^ but rather have 
the freedom of a variation of the moduli (pj under the condition cfp' = —r > 0. The 
two-dimensional noncollinear phase exists but has a marginal stability ^T]. AU other 
noncollinear phases even in the présence of a crystal anisotropy [v 7^ 0) either vanish 
or are unstable; for détails, see Réf. jTTj. This discussion demonstrates that the crystal 
anisotropy stabilizes the ordering along the main crystallographic directions, lowers the 
level of degeneracy of the ordered state related with the spontaneous breaking of the 
continuons symmetry and favors the appearance of noncollinear phases. 

The crystal field effects related to the unconventional superconducting order were estab- 
lished for the first time in Réf. ^Hj- In our considération of unconventional ferromag- 
netic superconductors in Sec. 4-7 we shall take advantage of thèse effects of the crystal 
anisotropy. In both cases f = and f 7^ the matrix (14) indicates an instability of 
three-dimensional phases (ail 0j 7^ 0) with an arbitrary ratios 0i/0j. As already men- 
tioned, for t> 7^ the phases of type 0i = 02 7^ 03 are also unstable whereas for t> = 0, 
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even the phase 0i = 02 = 03 > is unstable. 



4 M-triggered superconductivity 

We shall consider the Walker-Samokhin model |17] when only the M0i02— couphng be- 
tween the order parameters -0 and M is taken into account (7 > 0, 71 = 0) and the 
anisotropies (w = v = 0) are ignored. The uniform phases and the phase diagram in this 
case were investigated in Refs. [lH |13 HH] • 



Table 1: Phases and their existence and stabihty properties [6 = {62 — Oi), k = 0, ±1, ...]. 



N 


(f)j = M = 


always 


t > 0,r > 


FM 


= 0, = -t 


t < 


r > 0, r > re(t) 


SCI 


01 = M = 0, 02 = -r 


r < 


unstable 


SC2 


02 = = T:k, M = 


r < 


(t > 0)* 


SCS 


01 = 02 = M = 0, 0| = -r 


r < 


r <0,t>0 


COI 


01 = 02 = 0, 

02 = _r, m2 = -t 


r < 0, 
i < 


r < 
t < 


C02 


01 = 0, 02 = -r 

9 = 02 = TTk, M2 = -t 


r < 
t < 


unstable 


FS 


20f = 20| = 02 = -r + 7M 
03 = 0, 6» = 2TT{k - 1/4) 
7r = (72 - 2t)M - 2M^ 


7M > r 


3M2 > (-t + 72/2) 
M > 


FS* 


202 = 202 = 02 = -(r + 7M), 
03 = 0, = 2TT{k + 1/4) 
7r = (2t - 72)M + 2M3 


-7M > r 


3M2 > {-t + 72/2) 
M < 



Here we summarize the main results in order to make a clear comparison with the results 
presented in Sec. 5 and Sec. 6. Our main aim is the description of a trigger effect which consists 
of the appearance of a "compehed superconductivity" caused by the présence of ferromagnetic 
order (here, this is a standard uniform ferromagnetic order); see also Refs. |441 1451 EU] where 
this effect has been already established and briefiy discussed. As mentioned in the Introduction, 
a similar trigger effect is known in the physics of improper ferroelectrics. We shall set ^3 = 
and use the notation 9 = A0 = (02 — 9i). 

4.1 Phases 

The possible (stable, metastable and unstable) phases are given in Table 1 together with the 
respective existence and stability conditions. The normal or disordered phase, denoted in Table 
1 by N, always exists (for ail températures T > 0) and is stable for t > 0, r > 0. The 
superconducting phase denoted in Table 1 by SCI is unstable. The same is valid for the phase of 
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Figure 1: h — as a function of M for 7 = 1.2, and t — —0.2. The parameters r, t, 
and 7 are given by Eq. (11). 

coexistence of ferromagnetism and superconductivity denoted in Table 1 by C02. The N-phase, 
the ferromagnetic phase (FM), the superconducting phases (SCl-3) and two of the phases of 
coexistence (COl-3) are generic phases bccausc they appear also in the decoupled case (7 = 0). 
When the Mç!)i02-couphng is not présent, the phases SCl-3 are identical and represented by the 
order parameter (f) with components (f)j that participate on equal footing. The asterisk attached 
to the stabihty condition of the second superconductivity phase (SC2) indicates that our analysis 
is insufïicient to détermine whether this phase corresponds to a minimum of the free energy. 

It wih be shown that the phase SC2, two other purely superconducting phases and the coex- 
istence phase COI, have no chance to become stable for 7 7^ 0. This is so, because the phase 
of coexistence of superconductivity and ferromagnetism (FS in Table 1) that does not occur 
for 7 = is stable and has a lower free energy in their domain of stability. A second domain 
(M < 0) of the FS phase is denoted in Table 1 by FS*. Here we shall describe only the first 
domain FS. The domain FS* is considered in the same way. 

The cubic équation for magnetization of FS-phase (see Table 1) is shown in Fig. 1 for 7 = 1.2 
and t = —0.2. For any 7 > and t, the stable FS thermodynamic states are given by r(M) < 
fm = for M > Mm > 0, where corresponds to the maximum of the function r(M). 

The dependence of Mm{t) and Mo{t) = {-t + = ^/3Mmit) on t is drawn in Fig. 2 for 

7 = 1.2. Functions rm{t) = 4M^(i)/7 for t < 7^/2 (depicted by the line of circles in Fig. 3) and 



re{t)=7\t 



1/2 



(15) 



for i < define the border lines of stability and existence of FS. 
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Figure 2: The magnetization M versus i for 7 = 1.2: the dashed line represents Mq, the 
solid line represents M^q, and the dotted line corresponds to M^. 

4.2 Phase diagram 

We have outlined the domain in the (t, r) plane where the FS phase exists and is a minimum of 
the free energy. For r < the cubic équation for M (see Table 1) and the existence and stability 
conditions are satisfied for any M > provided t>'~f^. For t < 7^ the condition M > Mq have 
to be fulfilled, here the value Mq = {-t + 7^/2)^/2 Qf ^ is obtaincd from r(Mo) = 0. Thus 
for r = the N-phase is stable for t > 7^/2, and FS is stable for t < 7^/2. For r > 0, the 
requirement for the stability of FS leads to the inequalities 



where Mm = (Mo/\/3) and Mq should be the positive solution of the cubic équation of state 
from Table 1; > gives a maximum of the function r{M); see also Figs. 1 and 2. 

The further analysis defines the existence and stability domain of FS below the line AB denoted 
by circles (see Fig. 3). In Fig. 3 the curve of circles starts from the point A with coordinates 
(7^/2, 0) and touches two other (solid and dotted) curves at the point B with coordinates 
{ts = —7^/4, rs = 7^/2). Line of circles represents the function r{Mm) = rm{t) where 



Dotted line represents re(t) defined by Eq. (15). The inequality r < rm{t) is a condition for the 
stability of FS, whereas the inequality r < î'e(t) for (— t) > 7^/4 is a condition for the existence 
of FS as a solution of the respective équation of state. This existence condition for FS is obtained 
from 7M > r (see Table 1). 




(16) 




3/2 



(17) 
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Figure 3: The phase diagram in the plane {t, r) with two tricritical points (A and B) 
and a triple point C; 7 = 1.2. The parameters r ~ [T — Ts{P)] and t ~ [T — Tf{P)] are 
defined by Eq. (11). The domains of existence and stability of the phases N, FM and FS 
are shown. The line of circles represents the function rm{t) given by Eq. (17). The dotted 
line represents the function re(t) given by Eq. (15). On the left of point B, the same 
dotted curve corresponds to a FM-FS phase transition of second order. The equilibrium 
lines of N-FS and FM-FS phase transitions of first order are given by the solid lines AC 
and CB, respectively. 
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In the région on the left of the point B in Fig. 3, the FS phase satisfies the existence condition 
7M > r only below the dotted Une. In the domain confined between the hnes of circles and 
the dotted line on the left of the point B the stability condition for FS is satisfied but the 
existence condition is broken. The inequality r > re{t) is the stability condition of FM for 
< i-t) < 7^/4. For (-t) > 7^/4 the FM phase is stable for ah r > re(t). 

In the région confined by the line of circles AB, the dotted line for < (— t) < 7^/4, and the 
axis, the phases N, FS and FM have an overlap of stability domains. The same is valid for 
FS, the se phases and COI in the third quadrant of the plane {t, r). The comparison of the 
respective free énergies for r < shows that the stable phase is FS whereas the other phases are 
metastable within their domains of stability. 

The part of the t-axis given by r = and t > 7^/2 is a phase transition line of second order 
which describes the N-FS transition. The same transition for < t < 7^/2 is represented by 
the solid line AC which is the equilibrium transition line of a first order phase transition. The 
equilibrium transition curve is given by the function 

1 



r, 



eq\ 



Hère 



1 



„2 , i«^\l/2 



(0 = 4 37-(7' + 16t) ' M,,{t). (18) 



Me,{t) = ^ h' - 8i + 7 (7' + 16t)'^'j " " (19) 

is the equilibrium jump of the magnetization. The order of the N-FS transition changes at the 
tricritical point A. 

The domain above the solid line AC and below the line of circles for t > is the région of a 
possible overheating of FS. The domain of overcooling of the N-phase is confined by the solid 
line AC and the axes (t > 0, r > 0). At the triple point C with coordinates [0, î'eq(O) = 7^/4] 
the phases N, FM, and FS coexist. For t < the straight line 

7^ 

<q{t) = \ + \tl tB<t<0, (20) 

describes the extension of the equilibrium phase transition line of the N-FS first order transition 
to négative values of t. For t < ts the equilibrium phase transition FM-FS is of second order 
and is given by the dotted line on the left of the point B which is the second tricritical point in 
this phase diagram. Along the first order transition line r*q{t) given by Eq. ()2U() the equilibrium 
value of M is M(,q = 7/2, which implies an equilibrium order parameter jump at the FM-FS 
transition equal to (7/2 — y^)- the dotted line of the second order FM-FS transition the 
equilibrium value of M is equal to that of the FM phase {M^q = \/\t\)- The FM phase does not 
exist below Tg and this is a shortcoming of the model H12() with 71 = 0. 

The equilibrium FM-FS and N-FS phase transition lines in Fig. 3 can be expressed by the 
respective equilibrium phase transition températures T(,q defined by the équations = r{Teq), 
feq = r{Teq), r*^ = r{Teq), and with the help of the relation Meq = MiT^q). This limits the 
possible variations of parameters of the theory. For example, the critical température (Tgg = Te) 
of the FM-FS second order transition (7^/4 < —t) is obtained in the form Te = {Ts+^T^JAi/as)-, 
or, using M = (— a//6j)^/^, 

1/2 



1/2 



Tc = Ts + 

2 



(21) 
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Here Tj > T^, and T* = {'i-ïïjfaf/a'l 6j is a characteristic température of the model with 
= w = V = Q. A discussion of Eq. (21) is given in Sec. 5.3. 

The investigation of the conditions for the vahdity of Eq. 1)21(1 leads to the conclusion that the 
FM-FS continuons phase transition (at < —t) wih be possible only if the following condition 
is satisfied: 

Tf-Ts> = {ç + ^)T\ (22) 

where ç = bfa'^/Abga'j:. Therefore, the second order FM-FS transition should disappear for 
a sufïiciently large 7-coupling. Such a condition does not exist for the first order transitions 
FM-FS and N-FS. 

The inclusion of the gradient term (4) in the free energy (jj} should lead to a dépression of 
the equilibrium transition température. As the magnetization increases with the decrease of the 
température, the vortex state should occur at températures which are lower than the equilibrium 
température T^g of the Meissner state. For example, the critical température (Te) corresponding 
to the vortex phase of FS-type has been evaluated [17j to be lower than the critical température 
(PT|) : (Te — Te) = A-KfisM/as, where fis = \e\h/2mc is the Bohr magneton. For J S> fis, we 
have Te ^ Te- 

For r > 0, namely, for températures T > Tg the superconductivity is triggered by the magnetic 
order through the 7-coupling. The superconducting phase for T > T, is entirely in the (t, r) 
domain of the ferromagnetic phase. Therefore, the uniform supeconducting phase can occur for 
T > Tg only through a coexistence with the ferromagnetic order. 

The properties of the magnetic susceptibility and the spécifie heat near the phase transition 
lines shown in Fig. 3 have been investigated in Refs. |491l46j and here we shall not dwell on thèse 
topics. Note that the results [491 US] for the thermodynamic quantifies should be extended to 
include the physical effects considered in the next parts of this review. Such a considération 
requires a numerical analysis. 

In the next Sections we shall focus on the température range T > Tg which seems to be of 
main practical interest. We shall not dwell on the superconductivity in the fourth quadrant 
(t > 0, r < 0) of the {t, r) diagram where pure superconducting phases can occur for Systems 
with Tg > Tf, but this is not the case for UGe2, URhGe and ZrZn2. Also we shall not discuss 
the possible metastable phases in the third quadrant (t < 0, r < 0) of the {t, r) diagram. 

4.3 Note about a simplified theory 

The analysis in this Section can be done following an approximate scheme known from the theory 
of improper ferroelectrics; see, e.g., Réf. j55j . In this approximation the order parameter M is 
considered small enough which makes possible to ignore M^-term in the free energy. Then one 
easily obtains from the data for ES presented in Table 1 or by a direct calculation of the respective 
reduced free energy that the order parameters (/> and M of FS-phase are described by the simple 
equalities r = (7M — (fP') and M = (7/24)0^. For ferroelectrics working with oversimplified free 
energy gives a substantial departure of theory from experiment [HSl- The same approximation 
has been recently applied to ferromagnetic Bose-Einstein condensâtes [Hni EU • 

For ferromagnetic superconductors the domain of reliability of this approximation could be 
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the close vicinity of the ferromagnetic phase transition, i.e., for températures near the critical 
température Ty. This discussion can be worthwhile if only the primary order parameter also 
exists in the same narrow température domain (0 > 0). Therefore, the apphcation of the 
simphfied scheme can be useful in Systems, where >Tf. 

For Ts < Tf, the analysis can be simphfied if we suppose a relatively small value of the modulus 
(j) of the superconducting order parameter. This approximation should be valid in some narrow 
température domain near the line of second order phase transition from FM to FS. 

5 EfFect of symmetry conserving coupling 

Hère we shall include in our considération both linear and quadratic couplings of magnetization 
to the superconducting order parameter which means that both parameters 7 and 71 in free 
energy (12) are différent from zéro. In this way we shall investigate the effect of the symme- 
try conserving 71-term in the free energy on the thermodynamics of the System. When 7 is 
equal to zéro but 71 7^ the analysis is easy and the results are known from the theory of 
bicritical and tetracritical points jl3l 1531 1681 169j . For the problem of coexistence of conventional 
super conductivity and ferromagnetic order the analysis (7 = 0,71 7^ 0) was made in Réf. [à'2\ . 

Ai this stage we shall not take into account any anisotropy effects because we do not want to 
obscure the influence of quadratic interaction by considering too many parameters. For 7, 71 7^ 
and w = 0, V = the results again can be presented in an analytical form, only a small part of 
phase diagram should be calculated numerically. 

5.1 Phases 

The calculations show that for températures T > Tg, i.e., for r > 0, we have again three stable 
phases. Two of them are quite simple: the normal (A^-) phase with existence and stability 
domains shown in Table 1, and the FM phase with the existence condition t < as shown in 
Table 1, and a stability domain defined by the inequality r^^ < r. Here 

= lit + 7^/^, (23) 

and one can compare it with the respective expression (15) for 71 = 0. In this paragraph we shall 
retain the same notations as in Sec. 4, but with a superscript (1) in order to distinguish them 
from the case 71 = The third stable phase for r > is a more complex variant of the mixed 
phase FS and its domain FS*, discussed in Sec. 4. The symmetry of the FS phase coïncides with 
that found in j47j . 

We have to mention that for r < there are five pure superconducting (M = 0, <j) > 0) phases. 
Two of them, {(pi > 0, 02 = ç!)3 = 0) and {cpi = 0, 02 > 0, 03 > 0) are unstable. Two other 
phases, (0i > 0, 02 > 0, 03 = 0,02 = 9i + irk) and (0i > 0, 02 > 0, 03 > 0,02 = Oi + nk, 83 - 
arbitrary; A; = 0, ±1, ...) show a marginal stability for t > jir. 

Only one of the five pure superconducting phases, the phase SC3, given in Table 1, is stable. 
In case of 71 7^ the values of 0j and the existence domain of SC3 are the same as shown 
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in Table 1 for 71 = but the stability domain is différent and is given by t > 71 r. When 
the anisotropy effects are taken into account the phases exhibiting marginal stability within the 
présent approximation may become stable. Besides, three other mixed phases (M 7^ 0, > 0) 
exist for r < but one of them is metastable (for 7^ > l,i < 7ir, and r < 7it) and the other 
two are absolutely unstable. 

Here the thermodynamic behavior for r < is much more abundant in phases than for improper 
ferroelectrics with two component primary order parameter j53j . However, at this stage of 
expérimental needs about the properties of unconventional ferromagnetic superconductors the 
investigation of the phases for températures T < is not of primary interest and for this reason 
we shall focus our attention on the température domain r > 0. 

The FS phase for 71 / is described by the following équations: 

,/.i = (/>2 = ^, ^3 = 0, (24) 

(/>2 = (±7M-r-7iM2), (25) 
(1 - jj)M^ ± I771M2 _^f^_f__ M ± ^ = 0, (26) 



2 ' ' V 2 

and 

(^2-^1) = T| + 27rfc, (27) 

{k = 0, ±1, ...). The upper sign in Eqs. (24) - (27) corresponds to the FS domain where sin(02 — 
61) = —1 and the lower sign corresponds to the FS* domain with sin(02 — ^1) = 1- This is 
a generalization of the two-domain FS phase discussed in Sec. 3. The analysis of the stability 
matrix (14) for thèse phase domains shows that FS is stable for M > and FS* is stable for 
M < 0, just like our resuit in Sec. 4. As thèse domains belong to the same phase, namely, have 
the same free energy and are thermodynamically équivalent, we shall consider one of them, for 
example, FS. 



5.2 Phase stability and phase diagram 

In order to outline the (t, r) phase diagram we shall use the information given above for the 
other two phases which have their own domains of stability in the {t, r) plane: N and FM. The 
FS stability conditions when 71 7^ become 

27M - r - 71 > 0, (28) 

7M > 0, (29) 

3(1 - 7?)m2 + 3771M + t-jir- 7V2 > 0. (30) 

and we prefer to treat Eqs. (28) - (30) together with the existence condition 0^ > 0, with (p 
given by Eq. (25), with the help of the picture shown in Fig. 4. 
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The most direct approach to analyze the existence and stability of FS phase is to express r as 
of function of (M, t) from the équation of state (26), 



ril\t)= (31) 



eq 

(1 - ^f)M'' + ^771 + (t ■ 



Meq 



(7iMe, - 7/2) ■ ■ ■ 



and to substitute the above expression in the existence and stability conditions of FS-phase. It 
is obvions that there is a spécial value of M 

Msi = ^ (32) 
271 

that is a solution of Eq. (26) for any value of r and 

72 

for which this procédure cannot be applied and should be considered separately. Note, that 
M51 is given by the respective horizontal dashed line in Fig. 4. The analysis shows that in the 
interval ^ <t< 7^/2 the phase transition is again of first order; here 

- t (u) 

- 4(l + 7i)2- ^ ^ 

To find the equilibrium magnetization of first order phase transition, depicted by the thick line 
ACB in Fig. 4 we need the expression for equilibrium free energy of FS-phase. It is obtained 
from Eq. (12) by setting {w = 0,v = 0) and substituting r, as given by Eqs. (24), (25) and 
(31). The resuit is 

2 2 

+2[t{l - 7?) - ^]m2 - 277itM + t{t - ^)}, 

where M = Meq. 

For the phase transition from N to FS phase (0 < t < 7^/2), Meq is found by setting the FS free 
energy from the above expression equal to zéro, as we have by convention that the free energy 
of the normal phase is zéro. The value of Mel^ for positive t is obtained numerically and is 
illustrated by thick black curve AC in Fig. 4. When t^'* < t < the transition is betwcen FM 
and FS phases and we obtain Mg^^ from the équation fps = fpM = (^^^/2), where fpM is the 
free energy of FM phase. The equilibrium magnetization in the above i-interval is given by the 
formula 

2(1 + 71)' ^^^^ 

and is drawn by thick line CB in Fig. 4. The existence and stability analysis shows that 
for r > the equilibrium magnetization of the first order phase transition should satisfy the 
condition M^^ < M^l^ < M^^\ 
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Figure 4: The dependence M{t) as an illustration of stability analysis for 7 = 1.2, 71 = 0.8 
and w = 0. The parameters of the theory (r, t, 7, 71, -u;, . . . ) arc defined by Eq. (11). The 
horizontal dashed lines represent the quantities M51 given by Eq. (32) and Ms2 = 2M51. 
The line of circles AS1S2 describes the positive solution of Eq. (31). The thick line 
AC gives the equilibrium magnetization for t > 0. The thick line BC represents the 
equilibrium magnetization for t < as given by Eq. (35). The dotted curve is the 
smaller positive solution of the stabihty condition (30). The thin solid hne BS1S2 is the 
magnetization M — ^/^. The arrow indicates the triple point C. A and B are tricritical 
points of phase transition. The point Si corresponds to the maximum of the curve (23) 
for t <0, and the point 5*2 corresponds to ri^\t) = in Eq. (23). 

By Mq'^ we dénote the positive solution of A^^M^q) = and its t-dependence is drawn in Fig. 4 
by the curve with circles. Mm^ is the smaller positive root of stability condition (30) and also 
gives the maximum of the function rig^(M); see Eq. (31). The function Mm^ is dcpicted by 
the dotted curve AB in Fig. 4. Whcn t^i < t < the existence and stability conditions are 
fulfilled if < M < Msi, whcre is the magnetization of ferromagnetic phase and is 
drawn by a thin black line on the left of point B in Fig. (4). Here we have two possibilities: 
r > for < M < M^^^ and r < for M^^^ < M < M si- To the left of tsi and t > ts2, 
where 

the FS phase is stable and exists for Msi < M < ^Z— t. Here r will be positive when Mq^^ < 
M < and r < for M^^^ > M > Msi- When t <ts2, M < and r is always négative. 

On the basis of the existence and stability analysis we draw in Fig. 5 the {t, r)-phase diagram 
for concrète values of 7 and 71 . As we have mentioned above the order of phase transitions is 
the same as for 71 = 0, see Fig. 3, Sec. 4. The phase transition between the normal and FS 
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Figure 5: The phase diagram in the {t,r) plane for 7 = 1.2, 71 = 0.8 and w = 0. The 
parameters of the theory (r, t, 7, 71, w,...) arc dcfincd by Eq. (11). The domains of 
stabihty of the phases N, FM and FS are indicated. A and B arc tricritical points of 
phase transitions separating the dashed hnes (on the left of point B and on the right of 
point A) of second order phase transitions from the sohd hne ABC of first order phase 
transitions. The FS phase is stable in the whole domain of the {t, r) below the solid and 
dashed hnes. The vertical dashed hne coinciding with the r-axis above the triple point C 
indicates the N-FM phase transition of second order. 
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phases is of first order and goes along the equilibrium line AC in the interval {tA = 7^/2 and 
tc = 0). The function riç^(t) is given by Eq. (31) with M^l^ from Fig. 4. 

N, FM, and FS phases coexist at the triple point C with coordinates t = 0, and r^^ = 7^/4(71 + 
1). On the left of C for < t < the phase transition line of first order ri^*{t) is found by 
substituting in Eq. (31) the respective equilibrium magnetization, given by Eq. (35). In resuit 
we obtain 

This function is illustrated by the line BC in Fig. 5 that terminâtes at the tricritical point B 
with coordinates from Eq. (34), and 

(1) _ 7^(2 + 71) . . 

'^"4(1 + 71)^- ^^^^ 

To the left of the tricritical point B the second order phase transition curve is given by the 
relation (23) . Hcrc the magnetization is M = and the superconducting order parameter is 
equal to zéro {(/) = 0). This line intersects t-axis at ts2 and is well defined also for r < 0. The 
function r^\t) has a maximum at the point (t^i, 7^/471); hère M = Msi- When this point is 
approached the second derivative of the free energy with respect to M tends to infinity. The 
resuit for the curves r^q (t) of equilibrium phase transitions (N-FS and FM-FS) can be used to 
define the respective equilibrium phase transition températures Tps- 

We shall not discuss the région, t > 0, r < 0, because we have supposed from the very be- 
ginning that the transition température for the ferromagnetic ordering Tj is higher then the 
superconducting transition température T^, as is for the known unconventional ferromagnetic 
superconductors. But this case may become of substantial interest when, as one may expect, 
materials with Tf <Ts may be discovered experimentally. 



5.3 Discussion 



The shape of the equilibrium phase transition Unes corresponding to the phase transitions N-SC, 
N-FS, and FM-FS is similar to that of the more simple case 71 = and we shall not dwell on 
the variation of the size of the phase domains with the variations of the parameter 71 from zéro 
to values constrained by the condition 7^ < 1. Wc shall draw the attention to the important 
qualitative différence between the equilibrium phase transition lincs shown in Figs. 3 and 5. 
The second order phase transition line r^it), shown by the dotted line on the left of point B in 
Fig. 3, tends to large positive values of r for large négative values of t and remains in the second 
quadrant (t < 0, r > 0) of the plane {t, r) while the respective second order phase transition 
line ri^\t) in Fig. 5 crosses the t-axis at the point ts2 and is located in the third quadrant 
(t < 0, r < 0) for ail possible values t < ts2- This means that the ground state (at K) of 
Systems with 71 = will be always the FS phase while two types of ground states, FM and FS, 
can cxist for Systems with < 7^ < 1. The latter seems more rcalistic when we compare thcory 
and experiment, especially, in ferromagnetic compounds like UGe2, URhGe, and ZrZn2 where 
the présence of FM phase is observed at very low températures and relatively low pressure P. 

The final aim of the phase diagram investigation is the outline of the (T, P) diagram. Important 
conclusions about the shape of the {T,P) diagram can be made from the form of the (t,r) 
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diagram without an additional information about the values of the relevant material parameters 
{as, ttf, ...) and their dependence on the pressure P. One should know also the characteristic 
température Tg , which has a lower value than the experimentally observed HHl HH HHl OOl 1^ 
phase transition température {Tps ~ IK) to the coexistence FS-phase. A supposition about the 
dependence of the parameters as and a/ on the pressure P was made in Réf. [IZ]. Our results 
for Tf ^ Tg show that the phase transition température Tps varies with the variation of the 
System parameters (as,aj, ...) from values which are higher than the characteristic température 
Ts down to zéro température. This is seen from Fig. 5. 

In Systems where a pure superconducting phase is not observed for températures T ~ or 
T ~ Tps, we can set ~ in Eq. (21). Neglecting Ts in Eq. (21) and assuming that 
{T*/Tf) «; 1 we obtain that = Tps ~ (r*r/)V2. Note that the first (r*/r/)i/2.correction 
to this resuit has a négative sign which means that a suitable dependence of the characteristic 
température T* on the pressure P may be used in attempts to describe the expérimental shape 
of the FM-FS phase transition line in the (T, P) diagrams of UGe2 and ZrZn2; see, for example, 
Fig. 2 in Réf. Fig. 3 in Réf. |27], Fig. 4 in Réf. [3U| . The expérimental phase diagrams 
indicate that Tf[P) is a smooth monotonically decreasing function of the pressure P and Tf[P) 
tends to zéro when the pressure P exceeds some critical value Pc ~ 1 GPa. Postulating the 
respective expérimental shape of the function Tf{P) one may try to give a theoretical prédiction 
for the shape of the curve Tpg. 

The lack of expérimental data about important parameters of the theory forces us to make some 
suppositions about the behavior of the function T*{P). The phase transition température Tps 
will qualitatively follow the shape of Tf{P) provided the dependence T*{P) is very smooth. 
This is in accord with the expérimental shapes of thèse curves near the critical pressure Pc 
where both Tf and Tps are very small. The substantial différence between Tf and Tps at lower 
pressure (P < Pc) can be explained with the négative sign of the correction term to the leading 
dependence Tps{P) ~ [T* {P)Tf{P)]^/'^ mentioned above and a convenient supposition for the 
form of the function T*{P). 

Eq. (21) présents a rather simplified theoretical resuit for Te = Tps because the effect of M^|'(/'p 
coupling is not taken into account. But following the same ideas, used in our discussion of 
Eq. (21), a more reliable theoretical prédiction of the shape of FM-FS phase transition line can 
be given on the basis of Eq. (23). With the help of the experimentally found shape of Tf{P) and 
the définition of the parameters r and t by Eq. (11) we can substitute T = Tps{P) in Eq. (23). 
In doing this we have applied the following approximations, namely, that ~ for any pressure 
P, Tps{Pc) ~ Tf{Pc) ~ and for substantially lower pressure {P < Pc), Tf{P) > Tps{P). Then 
near the critical pressure Pc, we easily obtain the transition température Tps ~ 0, as should 
be. For substantially lower values of the pressure there exists an expérimental requirement 
(Tps — Ts) <C (Tf — Tps). Using the latter we establish the approximate formula 

{Tf-Tps)=j%^^h!af. (40) 

The same formula for (Tf — Tps) can be obtained from the parameter ts2{PFs) given by Eq. (37). 
The pressure dependence of the parameters included in this formula defines two qualitatively 
différent types of behavior of Tps{P) at relatively low pressures {P <C -Pc): (a) Tps{P) ^ 
below some (second) critical value of the pressure (P^ < Pc), and (b) finite Tps{P) up to P 0. 
Therefore, we can estimate the value of the pressure P^. < Pc in UGe2, where Tps{P'c) ~ 0. It 
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can be obtained from the équation 

Tj{P^) = i^Hf/j^af) (41) 

provided the pressure dependence of the respective material parameters is known. So, the above 
considération is consistent with the theoretical prédiction that the dashed hne in Fig. 5 crosses 
the axis r = and for this reason we have the opportunity to describe two ordered phases at 
low températures and broad variations of the pressure. Our theory allows also a description of 
the shape of the transition hne TpsiP) in ZrZn2 and URhGe, where the transition température 
Tps is finite at ambient pressure. To avoid a misunderstanding, let us note that the diagram in 
Fig. 5 is quite gênerai and the domain containing the point r = of the phase transition line 
for négative t may not be permitted in some ferromagnetic compounds. 

Up to now we have discussed expérimental curves of second order phase transitions. Our analysis 
gives the opportunity to describe also first order phase transition lines. Our investigation of the 
free energy (12) leads to the prédiction of triple (C) and tricritical points {A and B); see Figs. 
3 and 5. We shall not consider the possible application of thèse results to the phase diagrams of 
real substances, for which first order phase transitions and multicritical phenomena occur; see, 
e.g., Refs. [HSlIZni) where first order phase transitions and tricritical points have been observed. 
The explanation of the phase transition lines in Refs. [651 l7Uj requires further theoretical studies 
that can be done on the basis of a convenient extension of the free energy (12). For example, 
the investigation of vortex phases in Réf. jTO] requires to take into account the gradient terms 
(4). Another generalization should be done in order to explain the observation of two FM 
phases (nSUZOl- Note, that the experimentalists are not completely certain whether the FS phase 
is a uniform or a vortex phase, and this is a crucial point for the further investigations. But we 
find quite encouraging that our studies naturally lead to the prédiction of the same variety of 
phase transition lines and multicritial points that has been observed in récent experiments [651 

EOl. 



6 Anisotropy 

Our analysis demonstrates that when the anisotropy of Cooper pairs is taken into account, 
there will be no drastic changes in the shape the phase diagram for r > and the order of 
the respective phase transitions. Of course, there will be some changes in the size of the phase 
domains and the formulae for the thermodynamic quantifies. It is readily seen from Figs. 6 
and 7 that the température domain of first order phase transitions and the température domain 
of stability of FS above Tg essentially vary with the variations of the anisotropy parameter w. 
The parameter w will also insert changes in the values of the thermodynamic quantifies like the 
magnetic susceptibility and the entropy and spécifie heat jumps at the phase transition points. 

Besides, and this seems to be the main anisotropy effect, the w- and u-terms in the free energy 
lead to a stabilization of the order along the main crystal directions which, in other words, means 
that the degeneration of the possible ground states (FM, SC, and FS) is considerably reduced. 
This means also a smaller number of marginally stable states. 

The dimensionless anisotropy parameter w = Us/{hs + Us) can be either positive or négative 
depending on the sign of Ug- Obviously when Ug > 0, the parameter w will be positive too and 
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Figure 6: Phase diagram in the (i, r) plane for 7 = 1.2, 71 = 0.8, and w — 0.4. The 
meaning of hnes and points is the same as given in Fig. 5. 
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Figure 7: Phase diagram in the (t, r) plane for 7 = 1.2, 71 = 0.8, and w = —2. The 
straight dotted line for r < indicates an instability of the FS phase. The meaning of 
other lines and notations is the same as given in Fig. 5. 
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will be in the interval < w< 1 to ensure the positiveness of parameter b from Eq. (10). When 
u; < 0, the latter condition is obeyed if the original parameters of free energy (3) satisfy the 
inequahty —bs<Us<0. 

We should mention here that a new phase of coexistence of superconductivity and ferromag- 
netism occurs as a solution of Eqs. (13). It is defined in the following way: 



7? L 



71 U + 



2w 



— r 



(42) 



7? L 



71 r 



2w 



and 



2w sm{e2 -ei)= jM, cos(02 - ^i) 7^ 0. 



(43) 



(44) 



In the présent approximation the phase (42) - (44) is unstable, but this may be changed when 
the crystal anisotropy is taken into account. 



We shall write the équations for order parameters M and 
the changes when w 

2 ±7M - r - 71 



(l-w) 



<j of FS phase in order to illustrate 
> 0, (45) 



and 



{l-w- 7? )M^ 

3 o 
±-77iM2 ^ 



(46) 



t{l - w) 



71 

2 



71 r 



M± 



7r 



0, 



where the meaning of the upper and lower sign is the same as explained just below Eq. (27). 
The différence in the stability conditions is more pronounced and gives new effects that will be 
explained further, 

(2 - w)'yM - r - 71M2 



1 — w 



>0, 



(47) 



and 



7M -wr- wjiM^ > 0, 



3(1 - u; - 7^)M2 + 3771M + t{l -w)- 7^/2 - 7ir 

l-w 



> 0. 



(48) 
(49) 



The calculations of the phase diagram in {t, r) parameter space are done in the samc way as 
in case oî w = and show that for w > there is no qualitative change of the phase diagram. 
Quantitatively, the région of first order phase transition widens both with respect to t and r 
as illustrated in Fig. 6. On the contrary, when w < the first order phase transition région 
becomes more narrow but the condition (47) limits the stability of FS for r < 0. This is seen 
from Fig. 7 where FS is stable above the straight dotted line for r < and t < 0. So, purely 
superconducting (Meissner) phases occur also as ground states together with FS and FM phases. 
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7 Conclusion 



We investigated the M-trigger effect in unconventional ferromagnetic superconductors. This 
effect arises from the M^i^2-coupling term in the GL free energy and provokes the appearance 
of superconductivity in a domain of the system's phase diagram that is entirely occupied by the 
ferromagnetic phase. The coexistence of unconventional superconductivity and ferromagnetic 
order is possible for températures above and below the critical température Tg, that corresponds 
to the standard second-order phase transition from normal to Meissner phase - usual uniform 
superconductivity in a zéro external magnetic field which occurs outside the domain of existence 
of ferromagnetic order. Our investigation is mainly intended to clarify the thermodynamic 
behavior at températures Tg < T < Tf where the superconductivity cannot appear without the 
mechanism of M-triggering. We describe the possible ordered phases (FM and FS) in this most 
interesting température interval. 

The Cooper pair and crystal anisotropies are investigated and their main effects on the ther- 
modynamics of the triggered phase of coexistence is established. Of course, in discussions of 
concrète real materials the respective crystal symmetry should be considered. But the low sym- 
metry and low order (in both M and ip) 7-term in the free energy détermines the leading features 
of the coexistence phase and the dependence of essential thermodynamic properties on the type 
of crystal symmetry is not so considérable. 

Below the superconducting critical température a variety of pure superconducting and mixed 
phases of coexistence of superconductivity and ferromagnetism exists and the thermodynamic 
behavior at thèse relatively low températures is more complex than in improper ferroelectrics. 
The case Tf < Tg also needs a spécial investigation. 

Our results are referred to the possible uniform superconducting and ferromagnetic states. Vor- 
tex and other nonuniform phases need a separate study. 

The relation of the présent investigation to properties of real ferromagnetic compounds, such 
as UGe2, URhGe, and ZrZn2, has been discussed throughout the text. In thèse compounds the 
ferromagnetic critical température is much larger than the superconducting critical température 
(Tf ^ Tg) and that is why the M-triggering of the spin-triplet superconductivity is very strong. 
Moreover, the 71-term is important to stabilize the FM order up to the absolute zéro (0 K), as is 
in the known spin-triplet ferromagnetic superconductors. Ignoring 47] the symmetry conserving 
7i-term does not allow a proper description of the real substances of this type. More expéri- 
mental information about the values of the material parameters (ug, af, ...) included in the free 
energy (12) is required in order to outline the thermodynamic behavior and the phase diagram 
in terms of thermodynamic parameters T and P. In particular, a reliable knowledge about 
the dependence of the parameters ag and aj on the pressure P, the value of the characteristic 
température Tg and the ratio Ug/aj at zéro température are of primary interest. 
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